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We consider a biopolymer bundle consisting of filaments that are crosslinked together. The crosslinks are 
reversible: they can dynamically bind and unbind adjacent filament pairs as controlled by a binding enthalpy. 
The bundle is subjected to a bending deformation and the corresponding distribution of crosslinks is measured. 
For a bundle consisting of two filaments, upon increasing the bending amplitude, a first-order transition is 
observed. The transition is from a state where the filaments are tightly coupled by many bound crosslinks, to 
a state of nearly independent filaments with only a few bound crosslinks. For a bundle consisting of more than 
two filaments, a series of first-order transitions is observed. The transitions are connected with the formation 
of an interface between regions of low and high crosslink densities. Combining umbrella sampling Monte 
Carlo simulations with analytical calculations, we present a detailed picture of how the competition between 
crosslink shearing and filament stretching drives the transitions. We also find that, when the crosslinks 
become soft, collective behavior is not observed: the crosslinks then unbind one after the other leading to a 
smooth decrease of the average crosslink density. 
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I. INTRODUCTION 

The cytoskeleton is a complex meshwork made of long 
elastic filaments coupled together with the help of nu- 
merous, rather compact crosslinking proteins^. An im- 
portant aspect of cytoskeletal assemblies is their dynamic 
nature, which allows them to react to external stimuli 
and adapt their internal structure and mechanical prop- 
erties according to the needs of the cell. The reversible 
nature of crosslink binding is an important mechanism 
that underlies these dynamical processes. For exam- 
ple, living cells show complex rheological properties that 
range from fluidization to reinforcement under stress^ 3 -, 
and reversible bonds between cytoskeletal filaments have 
been proposed as key mechanisms in mediating between 
these contradicting behaviors-^. Similar effects are be- 
lieved to be important for the rheological properties of re- 
constituted F-actin networks^—, in particular at low fre- 
quencies that correspond to the lifetime of the crosslink- 
mediated bond2*i£. 

Another important class of cytoskeletal assemblies are 
filament bundles. Crosslinked F-actin bundles form pri- 
mary structural components of a broad range of cy- 
toskeletal structures including stereocilia, filopodia, mi- 
crovilli or the sperm acrosome. Type and properties 
of the crosslinking protein allow the cell to tailor the 
dimensions and mechanical properties of the bundles 
to suit specific biological functions. In particular, the 
mechanical properties of these bundles play key roles 
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in cellular functions ranging from locomotionii — to 
mechanotransductionM, and fertilization^. 

In-vitro experiments and modeling have emphasized 
the role of the crosslink stiffness in mediating bundle me- 
chanical properties^—. It is less clear, however, in how 
far the reversibility of the crosslinking bond may affect 
bundle mechanical or dynamical properties. On the one 
hand, one expects reversible bonds to influence the con- 
formational properties of the filaments. Examples for 
this dependency are the formation of kinks in the bundle 
contou r 19 ' 20 , or an unbundling transition as the binding 
affinity of the crosslinks is reduce d 21 ' 22 . Conversely, bun- 
dle conformational changes or the application of destabi- 
lizing forces 23 directly influence the binding state of the 
crosslinks. 

The aim of this article is to deepen our understand- 
ing of this complex interplay between reversible crosslink 
binding and bundle mechanical and dynamical proper- 
ties. We consider the nonlinear response of a reversibly 
crosslinked filament bundle to an imposed external force 
or deformation. In particular, we want to determine how 
the external driving is reflected in the internal degrees of 
freedom of the bundle, notably the binding state of the 
crosslinks. Combining simulations and theory we will 
show that, depending on the mechanical stiffness of the 
crosslinking agent, the fraction of bound crosslinks can 
display a sudden and discontinuous drop. This indicates 
a cooperative unbinding process that involves the crossing 
of a free energy barrier. Choosing the proper crosslink- 
ing protein, therefore, not only allows one to change the 
composite elastic properties of the bundle, but also the 
relevant time-scales: the latter can be tuned from the 
single crosslink binding rate to the (much longer) escape 
time over the free energy barrier (a short communica- 
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FIG. 1. (a) An end-grafted bundle of length L is brought from 
an initial unbend state (configuration 1) to a bend state (con- 
figuration 2). The bundle contour r±(z) is given by Eq. {TJ, 
where C sets the bending amplitude, (b) "Close-up" view of 
the bundle inner structure. In our model, the bundle consists 
of Nf filaments, each filament being a chain of N x beads con- 
nected by harmonic springs (horizontal bonds). The beads 
can slide along the bundle contour, as indicated by the dis- 
placement Uj, a . The filaments are joined to each other by 
crosslinks (vertical bonds) which may dynamically bind and 
unbind. 



tion of these results was recently presented by one of 
us in Ref. u4 ). We emphasize that related effects may 
be important in a variety of different biological contexts 
also. For example, cellular adhesion and locomotion are 
dependent on the formation of transient cell-to-matrix 
bonds 2 ^^. Simple theoretical models highlight the com- 
plex interplay between specific adhesion, mediated by the 
binding agent, unspecific interaction with the substrate, 
and cell-membrane elasticity 2 ^— . 

The outline of the paper is as follows. In Section [Til 
we define and motivate our bundle model. Next, in Sec- 
tion IIII1 we present our numerical results obtained in 
Monte Carlo simulations. In Section IIV1 the develop- 
ment of a theoretical model is described. The combined 
efforts of simulation and theory allow us to obtain a de- 
tailed understanding of the underlying physical mecha- 
nisms that affect the crosslink binding state. We end in 
Section [V] with a discussion of the implications and the 
experimental relevance of our results. 

II. MODEL 

We consider a bundle in a two-dimensional plane 
(Fig. Q^a)). The bundle has length L and, in the initial 
state, is oriented along the z-axis of a fixed laboratory 
frame (configuration 1). We now envision an experiment 
whereby the bundle is brought from the initial state to 
a "bent" state (configuration 2). In the bent state, the 
contour (shape) of the bundle is described by a trans- 
verse displacement r±(z), the functional form of which 
depends on the boundary conditions and the specific way 



of loading (e.g. bending or buckling). While the details 
of the loading are irrelevant for the subsequent analysis 
we choose 

9CT 2 

r ± {z) = — — (1 -cos(ttz/2£)) , (1) 

7T 

which mimics an experiment where an end-grafted bun- 
dle is deformed by a tip-load at its free end^i. Note 
that the choice of Eq. ((TJ) is not essential: qualitatively 
similar results are obtained with different bundle shapes 
also. The parameter C reflects the bundle curvature and 
will serve as a measure of the amplitude of the imposed 
bending deformation. 

In principle, crosslink reorganization may affect the 
local curvature and lead to the formation of kinks in 
the bundle contou r 19 ! 20 . In the following we are pri- 
marily interested in the effect of an imposed deforma- 
tion on the crosslink binding state, thus, bundle shape 
is assumed to be given and constant over the time- 
scale of interest. To bring the bundle from configuration 
1 — > 2 obviously requires a bending energy, Wbend oc 
Kf Jq (d 2 r±(z)/dz 2 ) 2 dz, with m the filament bending 
stiffness. However, as r±(z) is not allowed to change 
Wbcnd plays no role in what follows. 

The inner structure of the bundle is an array of 
i = 1, . . . , Nf parallel filaments that are spaced a dis- 
tance b apart (Fig. Hfb)). Each filament is a chain of 
a = X,...,N X beads, with harmonic springs (horizon- 
tal lines) joining nearest neighboring beads; the spring 
constant equals k s , and a denotes the equilibrium spring 
length. Beads can only slide along the contour, mak- 
ing their motions effectively one-dimensional. It there- 
fore suffices to assign a single number Ui iOL to each bead, 
denoting the relative displacement of that bead from 
its equilibrium position. The possibility of performing 
lateral motion transverse to the bundle axis is thereby 
neglected. At sufficiently low crosslink density the en- 
tropy stored in these bending degrees of freedom have 
been shown to drive an unbundling transition 21,22 . Here, 
we are primarily interested in highly crosslinked bundles 
away from the unbundling transition, such that the lat- 
eral degrees of freedom can be assumed to be frozen out. 

The filaments are joined to each other by crosslinks 
(vertical bonds in FigQJb)). A pair of beads can be 
crosslinked provided they are vertical nearest neighbors, 
i.e. a bead (i, a) can only be crosslinked to the two beads 
(i ± 1, a) and not to any other beads. The maximum 
number of crosslinks thus equals 

N max = (N t -l)N x , (2) 

but we emphasize that not all "allowed" vertical bonds 
are necessarily crosslinked, and so the actual number of 
crosslinks N will generally be lower. 

The two dominant contributions to the elastic energy 
of the bundle are an axial strain and a shear strain. The 
former is due to stretching of filaments and may be writ- 
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ten as 



N f JV X 
i=l a=2 



(3) 



with Ui, a the relative bead displacements, and k s the 
spring constant of the horizontal bonds dchned previ- 
ously. Note that k s is related to real material properties 
via k s = EA/a, where E is the filament Young modu- 
lus, A its cross-sectional area, and a the spacing between 
successive sites along the filament backbone. 

The resistance to shear deformations is mediated by 
the crosslinks: with their two heads crosslinks connect 
two neighboring filaments and provide a means of me- 
chanical coupling between them. While the form of the 
shear strain follows naturally from the basic definitions 
of continuum elasticity, it is nevertheless illustrative to 
discuss its physical basis. As a consequence of bundle 
deformation, encoded by r±(z) of Eq. ([TJ, filaments have 
to slip relative to each other, bringing the crosslinking 
sites out of registry and therefore leading to crosslink 
deformation. The slip is given by b6 a , where 



dr±(z) 



dz 



= aCN x sin(7ra/2iV x ), (4) 



z=aL/N x 



is the local tangent angle of the bundle at the site of 
the crosslink a. Bringing the crosslinking sites back into 
registry is possible if one of the filaments stretches out 
farther than its connected partner, Ui+i tC t = Ui^ a + b9 a , in 
order to compensate for the bending induced mismatch. 
The shear contribution to the elastic energy may thus be 
written as 



JV f -l N x 

f E E ; 



b6 n 



(5) 
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with k x the crosslink shear stiffness. In the above, n i a = 
(0, 1) are the crosslink occupation variables: ni M = 1 
means that between beads (i, a) and (i + l,a) a crosslink 
exists, while rii >a — means that no crosslink is present. 
A key ingredient of this work is the possibility of crosslink 
(un)binding: the crosslink occupation variables n^ a are 
therefore fluctuating quantities (in contrast to previous 
studie o 32 ' 33 where they were quenched). We emphasize 3 ^ 
that Eq. can also be derived by properly discretiz- 
ing the shear energy of an elastic continuum with shear 
modulus G = k x /a. 

The bundle model of the present study is thus defined 
by the Hamiltonian 



H„ + H y 



(6) 



i.e. the sum of the stretch and shear contributions. 
Hence, there is a competition between filament stretching 
and crosslink shearing which provides the fundamental 
physical mechanism that governs the phenomena to be 
described. In the sections to come, we will study Eq. © 



using mostly the grand canonical ensemble. That is, we 
fix the crosslink chemical potential /i, but the total num- 

is allowed to 



ber of crosslinks N 
fluctuate. 



EJVf-l y^- /v x 
i—l A^a — 1 



A. Units and conventions 

The key parameters in our model are the bending am- 
plitude C, the crosslink chemical potential /i, the number 
of filaments Nf, and the spring constants k s ,k x . We also 
introduce the crosslink density n = N/N max , with N the 
number of crosslinks between the filaments, and -/V max 
the maximum number possible (see Eq. ©). In what 
follows we choose /3k s = 100, but the ratio k x /k s will be 
varied (irrelevant factors of (3 — 1/ksT are thus absorbed 
in the spring constants). We also expect a dependence 
on the bundle length L — aN x , especially near phase 
transitions. The lattice constants a, b are set to unity. 
We consider an end-grafted bundle, corresponding to the 
boundary condition u^i = 0. All other displacements 
(ui. a> ±), as well as the bond occupation variables (rij iQ ), 
are fluctuating quantities. 



III. COMPUTER SIMULATION RESULTS 

The simulations are performed using grand canoni- 
cal Monte Carle 35 combined with an umbrella sampling 
scheme 3 ^ (see Appendix A). The key output is the (nor- 
malized) distribution P(N), defined as the probability to 
observe the bundle in a state with N crosslinks. The um- 
brella sampling scheme ensures that P(N) is measured 
over the entire range < N < A max , even in regions 
where P(N) is very small. As a consequence, in the vicin- 
ity of a first-order phase transition, our results are less 
susceptible to hysteresi s 37 ' 38 . 



A. The case JV f = 2 

We begin our simulations with a bundle consisting of 
Nf = 2 filaments. In Fig. [21 we plot the average crosslink 
density (n) — NP(N)/N niax , as function of the im- 
posed bending amplitude C for several values of k x /k s . 
In all cases, (n) decreases with C, showing that the 
crosslinks unbind upon bending. The striking feature is 
that, for k x /k s high enough, (n) drops extremely sharply 
at some special value of the bending amplitude. In this 
situation, the unbinding of crosslinks is a collective phe- 
nomenon, reminiscent of a first-order phase transition. 

We now specialize to k x /k s = 0.01, i.e. the largest 
value considered in Fig. O where the transition is dis- 
tinctly first-order. For a given bending amplitude C, we 
calculate the chemical potential fx* where the first-order 
transition occurs. To locate (J,*, is varied at fixed C un- 
til the "susceptibility" x = ((N 2 ) - (A) 2 )/A max reaches 
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FIG. 2. Behavior of the bundle upon bending. Plotted is 
the average crosslink density (n) versus bending amplitude 
C for fcx/fc s = 10" 5 , 10~ 4 , 10" 3 , 10" 2 (from top to bottom). 
For low values of k x /k s , (n) decreases smoothly with C; for 
higher values, a strong first-order transition is observed (data 
for N f =2,N X = 150, jit = 0). 



FIG. 4. lnP(n) at fi* of the susceptibility maximum (con- 
form Eq. (0). The left (right) peak corresponds to an un- 
bound (bound) bundle, while at intermediate densities the 
bundle is partially bound. Note the large free energy bar- 
rier AF separating the bound and unbound states (verti- 
cal arrow). The squares (a-e) indicate the densities n at 
which the crosslink profiles of Fig. [S] were measured (data 
for Nf = 2,7Vx = 300,fc x /fc s = 0.01, C = 0.195). 
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FIG. 3. Phase diagram of a bundle consisting of Nf = 2 fila- 
ments. Plotted is the chemical potential fi* of the first-order 
transition versus the bending amplitude (data for k x /k s = 
0.01, iV x = 300). 



a maximum, i.e. we numerically solve 

/i*(C) : x -> max - (7) 

The "phase diagram" of Fig. [3] shows pf versus C ob- 
tained in this way. This curve plays the role of a binodal: 
it separates the regime where the filaments are tightly 
bound by many crosslinks from the regime where they 
are only loosely coupled by few crosslinks. Note that the 
simulated binodal does not extend all the way C — >• 
but is "cut-off" at some threshold value C — Ct, which 
reflects the finite length of the bundle. Hence, Ct does 
not correspond to a critical point in the usual thermody- 
namic sensed. 

Next, we investigate how the crosslinks unbind at the 
transition (i.e. when p, = p,*). In Fig. 21 we show the 
logarithm of the distribution P{n) measured at p* (note 
that In P(n) may be regarded as minus the free energy 
of the bundle) . The distribution is distinctly bimodal, as 
is characteristic of a first-order transition^ .. In addition, 



we have checked that the barrier AF increases with the 
bundle lengths^, providing further confirmation that the 
transition is genuinely first-order—. 

To understand how the transition from the unbound to 
the bound state progresses, we associate features of the 
distribution P(n) to the spatial organization of crosslinks 
within the bundle. To this end, we introduce the crosslink 
density profile pi(a) = {rii^ a ) measured along the bundle 
contour a = 1,...,N X between "adjacent" filaments i 
and i + 1. Of course, for Nf = 2, there is only one such 
profile: pi(a). Some typical profiles are shown in Fig. [5j 
each obtained for a different value of the overall crosslink 
density n. In (a), we show pi(a) for n = 0.020, corre- 
sponding to the left peak in P(n) where the bundle is 
unbound. We observe that p\(a) > only in a small 
region near the grafted end, but rapidly decays to zero 
thereafter. Hence, there is an interface (domain wall) 
in the system, separating a region of high crosslink den- 
sity from one of low crosslink density. As n increases, 
the domain wall gradually shifts toward the center of 
the bundle (b,c) up to some threshold density n*. At 
n* , the domain wall jumps discontinuously toward the 
bundle end, yielding a constant crosslink density along 
the entire contour (d). The value of n* is given by the 
crosslink density where P(n) attains its minimum. Once 
the domain wall has "jumped" , increasing n further no 
longer affects the shape of the profile but merely raises 
the plateau value (e). 

The fact that the domain wall "jumps" at n* indicates 
another first-order transition. To make this explicit, we 
performed a number of canonical simulations (i.e. at fixed 
n), and measured the bundle free energy F(ai) as func- 
tion of the domain wall position ot\ (in simulations cut 
is set by the crosslink furthest away from the grafted 
end) . The result is shown in Fig. [5] for three values of 
the crosslink density n. In all cases, F(a\) reveals two 
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FIG. 5. Crosslink density profiles pi(a) for several values of 
the overall crosslink density n. The values of n in (a,e) coin- 
cide with the peak positions in P(n) of Fig. [4] and thus reflect 
the "coexisting" unbound and bound states. By increasing n, 
a domain wall gradually shifts toward the bundle center (b,c) 
up to a density n* where it "jumps" to the bundle end (d). 
For the parameters used here (Nf = 2, N* — 300, k x /k s = 
0.01, C = 0.195) we find n* ~ 0.36, which coincides with the 
crosslink density where P(n) attains its minimum. 
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FIG. 6. Bundle free energy F(ai) as function of the domain 
wall position ct\ for three values of the crosslink density n. As 
n increases, a first-order transition takes place at which the 
domain wall "jumps" from «i < N x /2 to a?i = 7V X (data for 
Nt = 2,N X = 300,fc x /fc s = 0.01, C= 0.195). 



minima: the minimum at a.\ < N x /2 (ai — N x ) cor- 
responds to the unbound (bound) state. Note that the 
overall shape of F(ai) rather resembles a cubic polyno- 
mial in oi, which is the standard form of the Landau free 
energy expansion to describe a first-order transition (ai 
being the order parameter, and n the temperature). For 
small n, the unbound state is stable (top curve). As n 
increases, a first-order transition takes place above which 
the bound state is stable (lower curve). Precisely at the 
transition, the minima have the same free energy, here at 
n w 0.31 (center curve). Note that this somewhat under- 
estimates n* ss 0.36 of Fig. [5] One reason for the discrep- 
ancy is the use of different ensembles (canonical versus 
grand-canonical) in systems of finite size. Another rea- 
son is that the (un)bound states remain meta-stable over 
a considerable range in n around the transition. To accu- 
rately obtain n* , one would need to perform an umbrella 
sampling simulation of the full two-dimensional proba- 
bility distribution P(N, a\). 



B. The case N f = 4 

More generally, for a bundle consisting of Nf filaments, 
we expect a sequence of Nf — 1 first-order transitions, one 
for each pair of adjacent filaments. To characterize these 
transitions, crosslink density profiles were measured for 
a bundle with Nf = 4. In this case, there are i — 1,2,3 
adjacent filament pairs, with corresponding density pro- 
files pi{a). The profiles are depicted in Fig. [7] for several 
values of the crosslink density n. For small n, the bundle 
is unbound: pi (a) is zero everywhere, except for a small 
region near the grafted end (n = 0.011). As n increases, 
the crosslinks preferentially bind the outer filament pairs 
(i = 1,3), while the center pair (i = 2) remains unbound 
(n = 0.039). Increasing n further, the initial symmetry 
between outer pairs gets broken: with equal probability, 
one of the outer pairs i = 1,3 is selected; the binding of 
crosslinks then continues for that pair only (n = 0.122) 
ultimately leading to the first transition of the sequence 
(ri = 0.211). After the first transition, we thus have a 
bundle where one of the outer filament pairs is bound, 
while the remaining two pairs are unbound. Next, the 
other outer pair begins to bind (n = 0.389) leading to 
the second transition (n = 0.444). We now have a bun- 
dle where both outer filament pairs are bound, and the 
"1 f> 3" symmetry is restored again. Not surprisingly, 
the third (and last) transition of the sequence involves 
the binding of the center filament pair. The mechanism 
is the same as before, featuring a domain wall (n = 0.800) 
that "jumps" at the transition (n — 0.944). 

Note that the binding transitions of the outer filament 
pairs occur relatively close to each other (the correspond- 
ing densities are ri\ ~ 0.2 and n-i ~ 0.4 ks 2n\, respec- 
tively). However, to induce the binding of the center pair, 
a significantly larger density ~ 0.9 is required. This 
becomes more pronounced in the (C, /i) phase diagram. 
To each transition / = 1, 2, 3 in the sequence corresponds 
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FIG. 7. Simulation evidence for the existence of a sequence of binding/unbinding transitions in a bundle consisting of Nf = 4 
filaments. Plotted are the crosslink density profiles pi(a) for several values of the crosslink density n. In each of the graphs, 
the solid curve corresponds to p2(a) of the center filament pair; dotted and dashed curves represent pi(a) and p3(a) of the 
outer pairs (data for iV x = 300, k x /k s = 0.01, C = 0.19). 
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FIG. 8. Phase diagram of a bundle consisting of Nf = 4 fila- 
ments (solid curves). Plotted is the chemical potential p*j of 
the / = 1, 2, 3 first-order phase transitions versus the bending 
amplitude; the deep lower curve corresponds to the binding of 
the center filament pair (data for k x /k s = 0.01, N x = 200). 
For completeness, the binodal for Nt = 2 is also shown (dot- 
ted curve). As in Fig. [3] the binodals do not extend to C — > 
due to finite bundle length. 



A. The case Nt = 2 without bending 

We first consider a two-filament bundle (Nf = 2) 
without external deformation (C — 0). Assume all 
crosslinks to be bound for the moment, ni j0l — 1, with 
a = 1, . . . , N x . The Hamiltonian of Eq. © then becomes 



r i Nx 

^bundle = — A l + o ^ 
a— 1 



|(A Q+1 -A Q ) 2 



k Y M 



with A a = U2, a — Uf. a and Y\ = k s /2. Strictly speaking, 
we also need to include the combination S Q = U2, a + 
u\. a . However, as the latter do not couple to the crosslink 
occupation variables, we need not consider them in our 
treatment. 

The idea is to iteratively integrate out the degrees of 
freedom A a , and to monitor the resulting "flow" of the 
coupling constant Y. Each time one of the A Q variables is 
integrated out, the Hamiltonian retains the above form, 
but with a renormalized coefficient Y given by the recur- 
sion relation 



a (local) maximum in the susceptibility, and so the tran- 
sition chemical potential \i\ can be evaluated via Eq. ([7]) 
as before. The resulting phase diagram now features 
three binodals, with those corresponding to the bind- 
ing of the outer filament occurring very close together 
(Fig. [8J. Note also that the first of these "outer" bind- 
ing transitions coincides with the binodal of the Nf = 2 
bundle. 



with the fixed-point 



2 k x +k s /2 + Y l 



( y/1 + 2k s /k x - l) . 



(8) 



(9) 



At the i-th iteration step, the partition function Z thus 
acquires a factor (k x + k s /2 + E^) -1 / 2 , such that one can 
write 



N x 

n 



-1/2 



(10) 



IV. THEORY 



We now present a theoretical description that captures 
most of the features observed in the Monte Carlo simu- 
lations. 



where we have used Eq. (j8|) and dropped an overall factor 

(k s /2)- N *'\ 

To see how the above generalizes to the case of open 
crosslinks let us assume that the crosslinks from sites 
a = j, . . . , j+lj — 1 are open. We will call this a "bubble" 
of length lj in the following. The associated variables 
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A Q can immediately be integrated over, with the effect 
of generating a term 



(kf /2)(A J+lj - Aj_x) , fcf = fc/fo + 1), 

in the renormalized Hamiltonian. Thus, to account for 
bubbles, we have to substitute the stretching stiffness k s 
with fcg ff in Eqs. © and ([TO)) , as well as to reinterpret N x 
as the number of bound sites: Abound = N. The resulting 
expression is an exact, albeit intractable, representation 
of the partition function. 

To make progress, we use a mean-field (MF) approx- 
imation, where we assume the crosslinks to be homoge- 
neously distributed along the bundle (the interface will 
be accounted for later). The actual bubble length thus 
gets replaced by its average value — > N/N 
where x denotes the fraction of bound crosslinks. Fur- 
thermore, we replace the renormalized coupling constant 
by its fixed-point value, Ti — > r^, which in our MF ap- 
proximation can be written as 

TooW = ^ (yi + 2xk s /k x - 1) . 

Within these approximations, the partition function of 
Eq. (fTO)) evaluates to 



t^n/n* 



-N/2 



(11) 



which can easily be evaluated numerically. The term 
Pn = ( jv ) represents the usual "entropy of mixing" and 
counts the number of crosslink configurations compati- 
ble with a given N = xN x . We will show in Appendix B 
how the averaged crosslink occupation (n) compares to 
our simulation results. We will furthermore show how 
one can improve the theory by explicitly incorporating 
bubbles using the "necklace model" of Fisher 43 i 44 . 



B. The case N{ = 2 with bending and interface 

Next, we incorporate a finite bending amplitude 6 a 
into the theory. In addition, as the simulations indicate 
the possibility of an interface between a region of high 
and low crosslink density, we must appropriately gener- 
alize the above MF approach to the crosslink occupation 
variables n\_ a . To this end, we assume the crosslinks to 
be homogeneously distributed in the region of high den- 
sity only. A sharp interface separates this region from 
one without crosslinks 



xN x ja.\ a < a>i, 







a > ati, 



(12) 



where ot\ 6 [0, Nx] is the (unknown) axial position of the 
interface. Below we will also use the normalized interface 
position 



As before, x denotes the fraction of bound crosslinks, 
which can be expressed in terms of the occupation vari- 
ables as x — J2 a n i,a/N x . 

We now make an "Ansatz" for the displacement de- 
grees of freedom Ui t<x . We assume that the quadratic 
terms in Eq. ([5]) are small whenever there are crosslinks 
that bind the two filaments together. That is, pro- 
vided ni^a — 1, the corresponding displacement A a = 
U2,a — ui,a °c a . In the region of low crosslink density 
we can assume that A Q = const. We thus obtain 



uosm(wa/2N x ) a < oi\, 
uq sin(7rai/2A r x ) a > or, 



(14) 



by requiring continuity at a = ct\ , and where also Eq. Q 
was used. Introducing these expressions into the Hamil- 
tonian of Eq. ([B]), and minimizing with respect to uq, we 
obtain the following "saddle-point" contribution to the 
effective free energy 



F sp (x,y) =N x Af c (y) 
with functions 



_ yfcjy) 
x f s {y) 



(15) 



2 f yL 

fc(y) = j^J cos 2 (its /2)Lds. 
f>(v) = jj" sm 2 (7rs/2L)d S . 



(16) 



and bundle length L = aN x . The relevant parameters 
are A ~ k s b 2 (Ca) 2 , which encodes the dependence on 
bending amplitude A oc C 2 , and x ~ (k s /k x )/N x rep- 
resenting the effects of the crosslink stiffness x oc l/k x . 
The point to note is that F sp (x,y) still depends on the 
crosslink occupation variable well as on the location 
of the interface ct\ (via the parameter y of Eq. p^|) ). 

In a previous report^! we presented the special case 
y = 1, i.e. without an interface being present. In this 
limit one obtains for the free energy 



F ap (x,y = l) = F ap (x)=N 



A 



1 + Xq/x 



(17) 



y = ai /N x e [0, 1] 



(13) 



This simple form, which is illustrated in Fig. [9j conveys 
an intuitive picture of how a finite bending amplitude 
may lead to a discontinuous reduction of crosslink oc- 
cupation x. When x <C xq the free energy essentially 
grows linearly, F sp (x) ~ xk Xl with the energy scale set 
by the crosslink stiffness k x . This indicates that each 
crosslink contributes a certain amount of deformation 
energy, while the filaments remain nearly undeformed. 
The few crosslinks present are just not strong enough 
to force the filaments into a deformed state. This sit- 
uation changes when i > io, where the free energy 
saturates at a value set by the filament stretching stiff- 
ness, F sp (x) ~ k s . Now there are enough crosslinks to 
stretch out the filaments, at the same time relieving their 
own deformation. Together with the binding enthalpy, 
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FIG. 9. Illustration of the mechanism that induces a discon- 
tinuous reduction of the crosslink occupation as function of 
bending amplitude and/or chemical potential. The top curve 
shows F sp (x) of Eq. (|17[l , with the dashed lines indicating the 
limiting behaviors for x —¥ and x — I 1. The binding enthalpy 
— xfi (lower line), which is to be added to F ap (x), "tilts" the 
free energy landscape, leading to a free energy (center curve) 
featuring a coexistence between two states (at high and low 
crosslink occupation x, respectively). As the control parame- 
ters A, fjj are varied, a discontinuous transition from one state 
to the other is thus observed. 



x=N/A/ x 

0.2 0.4 0.6 0.8 1 




0.2 0.4 0.6 0.8 1 
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FIG. 10. Contour plot of the total free energy F(x, y) given 
by Eq. (|18p . Note that only the region x < y is physically 
relevant as the total number of crosslinks N is always less 
than ai. The dots correspond to the state points (a-e) in 
Fig. [S] where the crosslink profiles were obtained in simula- 
tions. The inset shows the free energy F as function of the 
normalized interface position y for fixed x = 0.25,0.31,0.36 
(from top to bottom) , which compares well to the simulation 
result of Fig. [6] 



E = —x/j,, which leads to the usual tilting of the free en- 
ergy landscape, we obtain a total free energy that has two 
coexisting states, at high and low crosslink occupation. 
As the bending amplitude A or the chemical potential fi 
is varied, it is therefore possible to observe a discontinu- 
ous transition from one state to the other. 

The existence of an interface does not change fun- 
damentally this picture but adds a second reaction 
coordinate that the bundle can utilize in order to 
minimize its free energy during the unbinding pro- 
cess. Combining Eqs. (ITT|) and (|T5|) we obtain Z = 

E«*o£*=oe- mJV,ai) with 



-PF(N,ai) 



G) 



1 



-N/2 



x e 



(18) 

-f)(F sp +^N) 



which defines the total free energy F(N,a>i). The lat- 
ter is illustrated in Fig. [TU] using parameters that cor- 
respond to the discontinuous transition of Fig. [5] The 
figure strikingly shows the coexisting states at high and 
low crosslink density, and the formation of an interface 
upon passing over the intermediate saddle-point. The 
theory also reproduces the free energy barrier along lines 
of constant N (inset), in line with the canonical simula- 
tions of Fig. [5] The transition pathway followed in these 
simulations is indicated by the light (red) points. Af- 
ter passing the saddle-point, at x « 0.35, the interface 
"jumps" from the center of the bundle to the distant 
end. It is interesting to compare this "canonical" path- 
way with the general shape of the basin of attraction into 
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FIG. 11. (top) Average crosslink occupation (n) and (bottom) 
average interface location (ar) as function of bundle curva- 
ture C using Ax = 150, and for different crosslink stiffness 
k x /k s = 10~ 5 ,10" 4 ,4-10~ 4 ,10" 3 ,10~ 2 (from top to bottom). 
In both graphs, the lines are obtained without fit parameters 
from the theoretical model, Eq. (|18l) . while the symbols indi- 
cate simulation data. 
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the bound state. This seems to favor a pathway closer 
to the diagonal, with an associated continuous motion of 
the interface. 

From Eq. (fTHj) it is straightforward to calculate the av- 
erage crosslink density (n) , as well as the average location 
of the interface {ai). Both are compared to simulation 
results in Fig. llll and the agreement is remarkably good: 
as the bending amplitude increases, there is a discontin- 
uous jump in both the crosslink density, as well as in the 
interface position. 

C. The case N f > 2 

Let us now turn to bundles with more than two fil- 
aments. The simulations have indicated a sequence of 
Nf — 1 transitions, one transition for each adjacent fil- 
ament pair (Fig. [7]) . Upon increasing the bundle defor- 
mation C, the crosslinks in the central filament pair un- 
bind first. This unbinding transition leaves two smaller 
weakly coupled sub-bundles. The successive transitions 
then happen in the centers of these sub-bundles up until 
all filament pairs are unbound. 

As with each filament pair the number of variational 
parameters increases, a full theoretical treatment be- 
comes intractable. We therefore choose a semi-analytic 
treatment that explicitly accounts for crosslinks in the 
respective central layer only. The stretching degrees of 
freedom that do not belong to this layer are integrated 
out by assuming^ 

Ui+l,a ~ u i,a = A a , (19) 

independent of the layer index i. This indicates that fil- 
ament stretching increases approximately linearly with 
the filament index, i.e. with the distance from the cen- 
ter of the bundle (filaments farther out from the cen- 
ter "inherit" the stretching from all the filaments on 
the inside). Such a linear dependence constitutes a cen- 
tral assumption in classical continuum theories, such as 
Euler-Bcrnoulli or Timoshenko beam theories^. For 
sufficiently stiff crosslinks, and without considering the 
possibility of an interface (y — 1), this model can be 
mapped onto a two-filament bundle, conform Eq. (1171) . 
with effective iVf-dependent parameters, A(Nf) ~ 7V f 3 
and Xo(Nf) ~ Nf. We then use these Nf -dependent pa- 
rameters in the full free energy of Eq. (|18l) to calculate 
the average crosslink density for the given layer. 

Fig. H2] compares the result of this calculation (curves) 
with simulation data (symbols) for the case Nf = 4. The 
middle curve shows the average crosslink density (n) of 
the entire bundle versus bending amplitude C: the agree- 
ment with the simulations is quite remarkable. The lower 
curve shows the average crosslink density in the central 
layer, which unbinds at C « 0.05. Here, the theory 
slightly underestimates the simulation results, but the 
location of the transition is accurately reproduced. The 
upper curve shows (n) of the outer layer, which unbinds 
at a much larger bending amplitude C ~ 0.15. This curve 
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FIG. 12. The analogue of Fig. [2] but for Nf = 4. The key 
difference is that we now observe a sequence of unbinding 
transitions. Plotted is the average crosslink density (n) as 
function of the bending amplitude C for k x /k B = 10~ 3 , fi = 
0, iVx = 150 (curves show results obtained using the "Ansatz" 
Eq. (|19p ; symbols represent simulation data). The lower (up- 
per) curve shows the crosslink density in the central (outer) 
layer. The middle curve reflects the crosslink density of the 
entire bundle, which was obtained by adding the contributions 
from the individual layers. 



is nearly equivalent to that of a bundle with Nf = 2 fila- 
ments, in agreement with the binodal of Fig. [8] 



V. DISCUSSION 

We have studied the response of a reversibly 
crosslinked filament bundle to an imposed bundle defor- 
mation. The central quantity was the average crosslink 
density (n) and its dependence on the imposed curvature 
C of the bundle backbone. As compared to simple Lang- 
muir adsorption, (n) = 1/(1 + e /3M ) at chemical potential 
/it, one expects a decreasing crosslink occupation with in- 
creasing bundle deformation. The reason is that bundle 
deformation leads to a mismatch between the crosslink 
binding sites, and therefore to an elastic energy cost for 
binding. We found this basic mechanism to indeed be 
true, but the detailed phenomenology of crosslink un- 
binding turns out to be surprisingly rich and goes beyond 
a simple shift of the equilibrium state as would be char- 
acterized, for example, by an effective chemical potential. 

Our main result is the possibility of a cooperative and 
discontinuous reduction of crosslink occupation as the 
bundle deformation is increased. The reason for this be- 
havior is the competition between the energy scales of 
crosslink shear k x , and filament stretch k s , which is par- 
ticularly efficient when the crosslinks are stiff. An un- 
binding event will then affect the force balance in the 
filaments, with the potential of influencing the bundle 
state far away from the unbinding site. On the other 
hand, when fc x is small and the crosslinks soft, crosslinks 
unbind one after the other leading to a smooth decrease 
of the average crosslink occupation. 

We have characterized in detail the discontinuous un- 
binding transition and identified the existence of an inter- 
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face; the latter separates a region of high crosslink density 
from a region essentially free of crosslinks. The formation 
of an interface is a collective process in which crosslinks 
have to reorganize within the bundle and find new bind- 
ing sites. A similar effect has been discussed in the con- 
text of the twisting of helical filament bundles, where 
crosslinks organize into certain "binding zones''^. As one 
increases the bundle twist away from its preferred value, 
these binding zones become shorter and shorter, thus ne- 
cessitating collective reorganizations of many crosslinks 
simultaneously. As evidenced in Fig. [5J such a process 
implies the crossing of a free energy barrier. The associ- 
ated time-scale of escape over the barrier may be much 
longer compared to single crosslink (un)binding events. 

Such long time-scales may indeed be present in a re- 
cent experiment with F-actin bundles crosslinked by a- 
actinin^ 7 -. In the experiment a bundle was brought into 
a deformed configuration, where it was kept for either 
ten or 1000 seconds. After this waiting time the bundle 
was released and its relaxation was monitored. For the 
shorter waiting time the bundle showed the expected ex- 
ponential relaxation into the straight ground state. For 
the longer waiting time, however, the bundle did not re- 
lax back, but remained in a state with a considerable 
residual bending deformation. Apparently, upon bundle 
deformation, new binding sites become available that sta- 
bilize the bent shape by allowing the crosslinks to rebind 
in more favorable states that avoid crosslink straining. 
After releasing the bundle these crosslinks act to stabilize 
the bent contour, thus leading to a plastically deformed 
bundle, where the ground-state is no longer straight 48 . 
In line with our interpretation in terms of a free energy 
barrier, the apparent time-scale (the waiting time) neces- 
sary to observe bundle plasticity was considerably longer 
than the time required for (un)binding of the individual 
a-actinin linkers, which is on the order of seconds^. 

Strictly speaking, our model does not allow for plastic 
deformation as the crosslinks are assumed to always bind 
to the same, initial binding sites: a crosslink at site a only 
binds to the site a on the next filament. We explicitly ex- 
clude the binding of crosslinks between non-neighboring 
sites, e.g. between a and a ± 1. If the bending- induced 
mismatch, W a , between the original sites at a is large, 
those "new" sites may actually be more favorable in 
terms of crosslink energy. Binding to new sites may 
then help to "freeze-in" the applied bending deformation 
leading to a plastically deformed contour. Such a model 
would pose a challenging problem for a theoretical analy- 
sis and has to be left for future work. Instead, we propose 
a simple mapping that allows us to incorporate bundle 
plasticity into the "elastic" bundle model presented in 
this work. 

To this end, we assume the new binding sites to be 
optimal, in the sense that for the given imposed bundle 
contour no elastic energy cost is associated with the re- 
binding of a crosslink. Reshuffling a crosslink from its 
original position to a new site may then be conceived 
of as being an unbinding event at zero chemical poten- 



tial. For simplicity, we further assume that rebinding 
is fast enough, such that all crosslinks are either bound 
to original sites or to new sites. In this case, the num- 
ber of crosslinks bound to new sites, M, can be inferred 
without further calculation from the results presented in 
this work, (M) = 7V max — (A r ) AI= o. If we now release the 
bundle from its deformed state, filament elasticity will 
try to relax the bundle back to its original unbent state. 
The population of newly bound crosslinks, however, acts 
against this relaxation and stabilizes the bent contour. 
Within our assumption the additional contribution to the 
crosslink shear energy is 

Ny, 

tfx CW = ^£( A «+^«-0 2 > (20) 

X a—1 

where we assumed N{ = 2 for simplicity. Here, 0° is 
the tangent in the reference configuration at position a 
that was imposed during the waiting time, while 9 a cor- 
responds to the tangent acquired during the relaxation 
process. 

As can be seen, the magnitude of the stabilizing effect 
depends on M, which depends on the waiting time. If the 
crosslink is stiff enough, such that a free energy barrier 
is present, fast thermalization is prevented. In this case 
M ps on short time-scales and the bundle will behave 
elastically. On time-scales long compared to the escape 
time over the barrier, the number of crosslinks bound to 
new sites reaches its equilibrium value, M — (M), and 
the bundle is maximally plastic. In the experiment of 
Ref. H3 only two waiting times were accessible. It would 
be interesting to systematically change the experimental 
time-scale. One possibility could be to introduce a rate, 
at which bundle deformation is increased. In the context 
of protein unfolding, similar experiments^ have proved 
extremely useful to extract information on the underlying 
free energy landscape. 
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Appendix A: Monte Carlo method 

The bundle is represented by a two-dimensional (i = 
1, . . . , iVf) x (a — 1, ... , N x ) lattice. To each lattice site 
(i, a) a real number Ui i0l is attached, denoting the lo- 
cal axial displacement. In addition, occupation variables 
n-i.a = (0, 1) are attached to "vertical" nearest neighbor- 
ing pairs (i,a) and (i + I, a). We simulate in the grand 
canonical ensemble using a biased Hamiltonian 



W(N), 



(Al) 



with -Htmndie the "unbiased" Hamiltonian of Eq. (|6|) , and 
W(N) a weight function defined on the total number of 
crosslinks N. The purpose of W(N) is to sample all 
microstates < N < iV max with equal probability. To 
construct W(N), which is a priori unknown, we use suc- 
cessive umbrella sampling^. Once known, the sought-for 
probability distribution in the number of crosslinks fol- 
lows as P(N) oc exp(— )3W{N)). As Monte Carlo moves 
we use single bead displacements and crosslink bind- 
ing/unbinding moves, each attempted with equal prob- 
ability. In a displacement move, a lattice site (i,a) is 
selected randomly, and the current displacement w.; ;Q of 
that site is replaced by u • a = m_ a + S, with —0.1 < S < 
0.1 drawn uniformly randomly. The new displacement 
u[ a is accepted with probability Pdisp = rnin [l, e _/3Ai3 ] , 
where AH is the energy difference between initial and 
final state (since displacements do not change N, both 
Eq. © and Eq. (|A1[) can be used to compute the en- 
ergy difference). During a crosslink move, a vertical bond 
is selected randomly, and the corresponding occupation 
variable ni tOC is "flipped" (n^a — gets replaced by 1, and 
vice versa). The new state is accepted with probability 

-Pxlink = min [l )e -/3Aff blaa +/3 M AAT] ; with ^ the cross l in k 

chemical potential, AN = ±1 the change in the num- 
ber of crosslinks, and Aiibias the energy difference which 
must now be calculated using the biased Hamiltonian of 
Eq. HU). 



Appendix B: Comparison to simple mean-field 
theory and the necklace model 

The simple mean-field approximation developed in 
Eq. (fTTj) is capable of accurately describing the thermo- 
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FIG. 13. Average crosslink density (n) of an undeformed 
bundle (C = 0) consisting of Nf — 2 filaments as function of 
the chemical potential fi using k x /k B = 10~ 3 , iVx = 100. Note 
the logarithmic vertical scale! The thin curve is the solution 
of the one-crosslink model, as presented by Eq. (5) of Ref. [24], 
while the thick curves were calculated using Eq. (|11[) (lower 
curve) and the necklace model (middle curve). The latter 
indeed yields better agreement with the simulation results 
(dots). 



dynamic properties of the bundle without external de- 
formation, i.e. C = 0. One can slightly improve on this 
result by explicitly incorporating the bubbles along the 
lines of the classic "necklace model" 43 i 44 . 

The partition function for a bubble of length M is ob- 
tained from Eqs. © and (fTU|) by setting k x — \i — 



M 



Zo{M) = Y[(k s /2 + r< 



-1/2 



(1 



1=1 



M 



-1/2 



(Bl) 



where we defined a characteristic bubble size M c = 
k s /2T . The value of the coupling constant at the be- 
ginning of the bubble, To, can be taken equal to the 
fixed-point value in the neighboring bound segment, Tq — 
Too(^x)- The partition function of a bound segment of 
length TV is 



Zi(N) 



k > 



-1/2 



N 



(B2) 



The full partition function can be obtained from the gen- 
erating functions $ Q (z) = ^2^ = iZ N Z a {N) and the so- 
lution to the equation $o(z)<&i(z) — 1. From this the 
average crosslink density (n) is obtained in the usual way 
by differentiating with respect to /i. The quality of the 
different approximations is compared in Fig. 1131 



